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Abstract 

We study the low energy effective theory describing gravity with broken spatial diffeomorphism invari¬ 
ance. In the unitary gauge, the Goldstone bosons associated with broken diffeomorphisms are eaten and 
the graviton becomes a massive spin-2 particle with 5 well-behaved degrees of freedom. In this gauge, the 
most general theory is built with the lowest dimension operators invariant under only temporal diffeo¬ 
morphisms. Imposing the additional shift and SO(3) internal symmetries, we analyze the perturbations 
on a FRW background. At linear perturbation level, the observables of this theory are characterized 
by five parameters, including the usual cosmological parameters and one additional coupling constant 
for the symmetry-breaking scalars. In the de Sitter and Minkowski limit, the three Goldstone bosons 
are supermassive and can be integrated out, leaving two massive tensor modes as the only propagating 
degrees of freedom. We discuss several examples relevant to theories of massive gravity. 
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1 Introduction 


In general relativity, space-time diffeomorphism invariance is the local symmetry principle underlying 
gravitational interactions. One of most profound physical implications is the equivalence principle [1] [2] . 
However, on a specified space-time background, one or more of the diffeomorphisms are generally broken 
by gauge fixing, and the pattern of symmetry breaking constrains the low energy degrees of freedom 
and dynamics on that background. For instance, our expanding universe can be considered as a temporal 
diffeomorphism breaking system, because the future always looks different from the past. Theories of grav¬ 
ity with temporal diffeomorphism breaking have been extensively studied in the literature, e.g. k-essence 
[3], the effective field theory of inflation [4] [5], ghost condensation [6], Horava gravity [7], generalized 
Horndeski theories [8] [9] [10] [11] and so on. 

Spatial diffeomorphism breaking is also important for the description of our universe: an everyday 
example is the description of low energy excitations of solids (phonons), which can be derived as the theory 
of broken spatial diffeomorphism invariance in which the phonons are the Goldstone bosons [12] [13]. With 
the addition of a 17(1) symmetry to conserve particle number, the theory of broken spatial diffeomorphisms 
describes “superfluid solids” (“supersolids”) [14]. In these systems, spatial diffeomorphism invariance is a 
hidden symmetry that is evidenced by the constrained form of the Goldstone bosons’ interactions. At solar 
system and cosmological scales, spatial diffeomorphism invariance is a relevant symmetry in that these 
systems are accurately described by general relativity. However, the unexplained origins of inflation, the 
end of inflation and the late time accelerated expansion keep open the possibility that general relativity 
is modified at the largest and smallest length scales. It is therefore interesting to ask how broken spatial 
diffeomorphisms impact cosmological dynamics. 

In this work, we develop the effective theory for the long-wavelength (k/a ~ H the Hubble constant) 
degrees of freedom in the presence of broken spatial diffeomorphisms. As in the condensed matter ex¬ 
amples, spatial diffeomorphism invariance can be broken by non-gravitational interactions. Field theory 
provides a mechanism in the form of soliton field configurations, such as the hedgehog solution 

r = f(r)—, a = 1,2,3, (1.1) 

r 

which describes a monopole in an SU(2) gauge theory that is spontaneously broken down to 17(1). Here a 
is the internal index when it is written as the superscript of scalar fields and is the spatial index when the 
superscript of coordinates. Taking into account gravity, this configuration of the cj) a fields breaks spatial 
diffeomorphisms, and in this case, translation and rotation symmetry are also broken to subgroups by 
fixing a preferred origin of the monopole. This background field configuration has been implemented to 
produce an inflationary phase in a model known as “topological inflation”, given that the size of monopole 
is greater than the Hubble radius in the early universe [15] [16]. The field configuration Eq. (1.1) is not 
the unique way to break spatial diffeomorphisms, and we will consider a more minimal way that preserves 
the translation and rotation symmetries. 

The low energy description of broken spatial diffeomorphisms exhibits three Goldstone bosons, scalar 
fields 4> a which physically can be thought of as measuring spatial position. In unitary gauge, these “ruler 
fields” are identified with the coordinates, 

cj) a = x a , a = 1,2,3. (1.2) 
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Translation and rotation invariance are preserved by implementing a shift symmetry (j) a —> (j) a + c a for 
constants c a and an SO(3 ) internal symmetry in the triplet <j) a . The scalars </>“ select a frame of reference, 
a background against which to measure perturbations. To restore the Goldstone bosons as dynamical 
degrees of freedom, we add a fluctuating component to the field 

(,b a ^x a + -K a (1.3) 

with 7 T a transforming under spatial diffeomorphisms opposite to the spatial coordinates and thus furnishing 
a nonlinear realization of the symmetry (known as the Stuckelberg trick). To see how this describes a 
solid, think of the scalar functions cj> a (x ) as locating each volume element or lattice site in space. In the 
long-wavelength limit A 3> lattice spacing, inhomogeneity at the sites is smoothed over, and fluctuations 
of the cj) a correspond to fluctuations of the site locations, i.e. phonons [12, 13]. 

Broken spatial diffeomorphism invariance is interesting in the context of gravitational theory, because it 
generates a mass for the graviton. This is easy to understand seeing that the presence of a fixed frame (one 
may think of a lattice) admits the propagation of additional compressional and rotational modes, which 
are the longitudinal modes of the graviton. The structure of the broken spatial diffeomorphism theory 
thus helps understand how to construct a general self-consistent theory of massive gravity. Indeed, it is a 
basic question in classical field theory whether an analog of Higgs mechanism exists that can give gravitons 
a small but non-vanishing mass. Experimentally, we do not know how gravity behaves at distances longer 
than ~ 1 Gpc, and the extremely tiny energy-scale associated with the cosmic acceleration[17, 18] hints 
that gravity might need to be modified at such large scale. 

The theoretical and observational consistency of massive gravity has been a longstanding problem. In 
the pioneering attempt in 1939 by Fierz and Pauli [19], the simplest extension of GR with a linear mass 
term suffers from the van Dam-Veltman-Zakharov discontinuity [20] [21], giving rise to different predictions 
for the classical tests in the vanishing mass limit. This problem can be alleviated by introducing nonlinear 
terms [22]. However, in 1972, Boulware and Deser pointed out that a ghost generally reappears at the 
nonlinear level, which spoils the stability of the theory [23]. Inspired by effective field theory in the 
decoupling limit [24], people have learned that in principle the Boulware-Deser ghost can be eliminated 
by construction [25] [26]. This theory is now dubbed dRGT gravity. When we apply dRGT gravity 
to cosmology, a self-accelerating solution is found for the open FRW universe [27]. However, follow-up 
cosmological perturbation analysis found a new ghost instability among the 5 gravitational degrees of 
freedom [28] [29] [30] [31] [32], and further dRGT gravity might suffer from acausality problems [33] [34]. 
The dRGT ghost instability can be eliminated at the expense of introducing a new degree of freedom 
[35] [36]. In this context, it is interesting to search for a simpler and self-consistent massive gravity theory, 
as an alternative to the Fierz-Pauli family of theories. 

An alternative way to realize a massive gravity theory is to break the Lorentz symmetry of vacuum 
configuration, in addition to the space-time diffeomorphisms. A broad class of Lorentz-symmetry breaking 
massive gravity theories have been discussed in Refs. [37] [38] [39] [40] [41], Among these theories, a simple 
example is the spatial condensation scenario Eq. (1.2); the non-vanishing spatial gradient breaks 3 spatial 
diffeomorphisms, while temporal diffeomorphism, translational and rotational invariance are preserved 
[42] [43]. Previous analyses focused on linear theory in the decoupling limit. As we will see below, the 
theory becomes degenerate in the Minkowski space time. On FRW backgrounds, there are exactly 5 
degrees of freedom in the theory. In the unitary gauge, the graviton eats the Goldstone excitations ir a in 
Eq. (1.3) and becomes a massive spin-2 particle, with 5 massive modes in the spectrum. 
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The resulting theory has several interesting applications. For instance, the graviton mass removes the 
IR divergence in graviton scattering [43], and leaves an interesting imprint on CMB primordial tensor 
spectrum [44] . A viable massive gravity theory also provides the basis for holographic study of dissipative 
systems [45-47] . Several other gravitational phenomena associated with broken spatial diffeomorphisms 
have also been discussed in the literature [48] [49] [50] [51] [52] [53] [54] [55] [56], without relating them ex¬ 
plicitly to the massive gravity aspect of the theory. For example, by tuning the form of higher order 
interactions, Ref. [50] builds a model of inflation, calling it “solid inflation”, in which they calculate 
the two- and three-point correlation of primordial perturbation. Our analysis helps understand why the 
sound speeds of scalar and vector modes are related in such a theory. 

In this paper, we study the general gravitational action for broken spatial diffeomorphisms by con¬ 
structing the appropriate low energy effective held theory. The effective held theory approach describes 
a system through the lowest dimension operators compatible with the underlying symmetries. Usually, 
when we study a gravitational system, we hrst write down a general covariant action, and space-time 
diffeomorphisms are broken “spontaneously” after solving the equation of motion. However, in this pa¬ 
per to learn more of the structure of the theory and resulting character of massive gravity, we start by 
writing down the most general gravitational action compatible with spatial diffeomorphisms breaking in 
the unitary gauge. We then recover general covariance by performing a change of spatial coordinates 
x l —> x l + £* and promote the parameters £* to Goldstone bosons which transform opposite to the spatial 
coordinates n a —> tt 0, — £“ under spatial diffs. 

This paper is organised as follows: In section 2, before constructing the specific effective theory, we 
discuss the general set of terms allowable in unitary gauge. Because the unitary gauge action, in its 
initial background-independent form, does not make explicit the dynamical degrees of freedom, we must 
carefully select the terms so as to preserve the 5 desired degrees. Then in section 3, we specify to the FRW 
background, discuss the physical scales of interest, including different requirements during inflation and 
late-time and requirements for the perturbativity of the theory. Restricting to SO(3 ) rotational symmetry 
and shift symmetry, we determine the effective action in the FRW universe, and analyse all scalar, vector, 
and tensor degrees of freedom. In section V, we present several examples of the applications of our 
formalism. Conclusion and discussion will be in the final section VI. In this paper we use the (—, +, +, +) 
convention in the space time metric. 

2 Generic action in unitary gauge 

To help show how a theory of broken spatial diffeomorphisms is a theory of massive gravity, we first 
discuss constructing the Lagrangian in the unitary gauge, in which we only have metric degrees of freedom. 
When we analyze perturbations, we identify which metric components become dynamical, corresponding 
to the longitudinal polarizations of the graviton, and thus in a given allowed operator we can track the 
real degrees of freedom. This is important because much previous study of modified gravity theories has 
shown that many forms of higher derivative operators lead to new degrees of freedom. To preserve exactly 
5 dynamical degrees of freedom (2 graviton polarizations + 3 goldstone bosons), the set of operators must 
be additionally constrained. As these constraints apply to the construction in any background, it is worth 
investigating allowable terms in a “generic” action in unitary gauge, before specifying the background 
solution (and with it the relevant scales). This facilitates construction in other backgrounds. 
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In unitary gauge, the action is only invariant under the temporal diffeomorphisms. There is a preferred 
spatial frame generated by the space like gradient of scalar fields, d^ a d v (j) b hab > 0, where 0 a (t,x) 
generally is a function of space and time, and h a b is the internal metric of scalar fields’ configuration. In 
the unitary gauge, the spatial frame x a is chosen to coincide with c/) a , 

(r) = x a , a= 1,2,3. (2.1) 


They transform as the scalars under the residual diffeomorphisms, so that the additional degrees of 
freedom are in the space-time metric. 

Going systematically through the geometric objects, we have: 

1. Terms that are invariant under all diffeomorphisms. These include polynomials of the Riemann 
tensor R^, UK \ and its covariant derivatives contracted to give a scalar. However, many such terms 
introduce additional unwanted degrees of freedom and/or break temporal diffeomorphisms. For 
instance, by doing a conformal transformation, 1Z 2 is equivalent to Einstein gravity plus a scalar 
field with non-trivial potential. To remain within the effective theory and its degrees of freedom, 
we need only the linear term in Ricci scalar 7Z to the order considered. 


2. Any scalar function of coordinates x a , as well as their covariant derivatives. In the unitary gauge, 
V^x a = <5“. Higher than second order derivatives generally give rise to extra modes and the classical 
Ostrogradski ghost instability. However, provided we have a viable perturbative expansion, in which 
higher dimensioned operators including higher derivative operators are supporessed by a high scale 
A, the typical mass of the Ostrogradski ghost is at or above the cut-off scale of our effective field 
theory. For this reason, the would-be ghosty modes are non-dynamical and can be integrated out 
at the low energy scale [5] . Consider for instance the scalar field theory with higher order derivative 
terms 


£ 


d^d^ + A- 2 (n^ 2 , 


for which the propagator reads 


A(k) 


¥ + fc4 


A 7 


1 1 

¥ ~ ¥ + A 2 ’ 


( 2 . 2 ) 


(2.3) 


with two propagating degrees. The second appears to have the wrong sign propagator, which is 
the possible ghosty mode. However, the pole is at ¥ = —A 2 , which is around the cutoff scale of 
our theory. For momenta in the domain of the effective theory k <C A, this degree of freedom is 
supermassive and can be integrated out. 


3. We can leave free the upper indices i in every tensor. For instance we can use g 1 ^, 7Z 12 and 7 Z^ kl . 
However, we notice that 1Z IJ can be rewritten into higher order covariant derivatives of x a by partial 
integration, 


U aa = IZ^d^d^ = {Ux a f - (V^V^x 0 ) (\7 ti \7 u x a ) + total derivatives, (2.4) 

so is 7 Z^ kl term. Thus VJi and 7Z l i kl terms are irrelevant at low energy scale. 
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We conclude that the most generic Lagrangian in the unitary gauge is given by 

S = J d A x^gF , (2.5) 

where all the free indices inside of the function F must be i’ s. The construction so far is general to 
the extent of a gravitational action respecting time diffeomorphism invariance and excluding degrees of 
freedom additional to the broken spatial diffeomorphisms. Additional symmetries must be respected when 
the theory is considered on a particular background, and in the next section we will discuss the restriction. 

3 Expanding around a FRW background 
3.1 Scales, power counting and constraints 

In the present work, we are primarily interested in cosmological phenomenology of the broken spatial diff 
theory, and hence we are looking at wavelengths of the Hubble scale. To be consistent with the high 
degree of spatial isotropy and homogeneity at cosmological lengths, the effective scalars vr a have an SO(3) 
internal symmetry and a shift symmetry 7r a TT a + c a . The shift symmetry requires the Goldstones 
have only derivative couplings, and consequently the effective theory is essentially an expansion in k/M 
where M is a scale appearing in association with the higher dimensioned operators. For the theory to 
be perturbative in that lower order (lower derivative) operators are more relevant than higher orders, 
we should have H <C A the breakdown momentum scale, which is approximately the momentum scale 
at which 7T7 t scattering violates unitarity. A is related to the Ms appearing in various higher derivative 
operators and thus the condition ff < A yields constraints on the sizes of the higher order operators. 
Conversely, we constrain M considering (lack of) evidence for higher order operators. Satisfying these 
constraints, power counting derivatives k/M yields a consistent theory because the breakdown scale is 
parametrically above the dynamical range of the theory. 

Note that the nontrivial background, the FRW spacetime, introduces an important scale in the dy¬ 
namics. H and H enter by determining the characteristic scale of the background (j) a fields, and as a 
consequence their fluctuation components n a . As we will show in the subsection 3.3, in the H —> 0 limit, 
the theory becomes strongly coupled and encounters some of the issues well-known to massive gravity on 
Minkowski backgrounds, perhaps pointing to a deeper physics reason for these issues. On the other hand, 
it has been previously proven that the present theory is continuous with GR in the limit M —> 0, with the 
Goldstones decoupling, becoming nondynamical and restoring the full diffeomorphism symmetry [43]. 

Considering the early universe, the kinetic energy of a scalar field with k H is ~ y/2eM p \H where 
e = —H/H 2 is the slow roll parameter, meaning its change in amplitude over a Hubble time H~ 1 is 
~ y/2eM p \. Since the triplet of dynamical scalars ir a can at most be responsible for inflation, the mass 
scale suppressing higher order terms must be larger M > \/2eM p i [5]. As we will show, achieving a 
self-accelerating cosmology with a minimal form of this theory requires fine-tuning the parameters in 
such a way that the theory becomes strongly-coupled and loses its meaning. We can consider instead 
inflation by independent dynamics, with small ir a fluctuations on top. In this case, the scale of the vr a 
must be smaller than the inflation energy scale, but still large enough to validate an expansion in k/M ; 
this scenario is quite natural in the context of GUT-scale inflation and the possibility of topological 
defects arising from breaking the GUT symmetry group. In fact, the energy density of the scalar kinetic 
term dilutes as d^ir a d u ir a ~ a~ 2 the same as a network of topological defects [58] [59], making such 
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an identification tempting. Since the 7r a are small perturbations to dominant inflationary background, 
evidence of this scenario must be sought in the gravitational wave signal from inflation, via the (effective) 
mass for gravitational waves [44]. 

At late times z < a few, H is much smaller, and the mass scale of the spatial diff breaking mechanism 
certain to be higher, in particular if it is related to any standard model particle or astrophysics. Con¬ 
sidering the present universe, we can constrain the scales via the resulting mass and sound speed. For 
non-vanishing graviton mass, the orbits of binary system will decay at a slightly faster rate than pre¬ 
dicted by GR, due to the additional energy loss from in the emission of massive gravitational waves. The 
decay rate difference can be roughly estimated as <5 ~ (~f) , where u is the frequency of gravitational 
waves, which is identical to the inverse of orbital period of binary, and m g is the graviton mass nowa¬ 
days. Consider the Hulse-Taylor binary pulsar, PSR B 1913+16, for which the observed orbital decay 
attributed to gravitational wave emission agrees with the predictions of GR to 0.3%, 5 < 0.3% requires 
that nig < lCT 20 eE 10 12 Hq [60]. Non-vanishing graviton mass also changes the propagating speed of 
gravitational waves, which leads to the upper bound on the graviton mass of m g < 10 _23 eV"[61]. The 
detection of gravitational waves in advanced LIGO could then bound the graviton mass potentially all the 
way down to m g < lCP 29 eE [61] [62] [63]. See [64] for a recent review on the theoretical and experimental 
aspects of massive gravity. 


3.2 Construction of the Action 

For a given FRW, the enhanced symmetry reduces the set of allowed operators from those discussed in 
the previous secion. In addition to preserving the unbroken temporal diffeomorphism invariance, we must 
preserve homogeneity (translation invariance) and isotropy. Terms like x a without covariant derivative 
operators acting on them are not allowed, because they break homogeneity and isotropy. Instead, x a 
terms must appear in the action with derivatives, 

f gv = d g x a d v x b 5 ab , (3.1) 


and any pair of <9^x a should be contracted with internal metric 5 ab to maintain the 50(3) spatial rotation 
symmetry. Each d lx x a d v x b 5 ab is contracted with the metric g lw and thus gives rise to terms proportional 
to the trace of spatial metric g n and cross terms like g 13 g 13 in the action. At linear perturbation level, 
the cross terms decompose into the trace sector, which is a function of g n , plus the traceless sector which 
only appears at quadratic order in the action for perturbations. 

As for higher derivative terms, they are “less relevant” at low energy scale, as we discussed in the 
section 2. Thus we only focus on lowest dimensional operators for the time being. The higher order 
derivative terms will be informatively discussed in section 3.5. 

We can now write down our most generic action in the unitary gauge as follows, 
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l -Mln + K + cY J g ii + - 

i 


(3.2) 


where the dots stand for terms which are of at least quadratic order in the fluctuations. The terms g n is 
responsible for the “spatial diffeomorphisms breaking”. Notice that the above three terms contain linear 
perturbations around FRW background. Therefore, the coefficients A and c will be fixed by requiring 
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that all tadpole terms cancel around the given FRW solution. The differences between models will be 
encoded into higher order terms. 

To fix the coefficients A and c of linear terms, we insert the FRW ansatz for the background space 
time metric, 

ds 2 = — dt 2 + a 2 dx 2 . (3.3) 

The Einstein equations read 

3M 2 H 2 = -34 - A, 
y a- 

M?H = (3.4) 

Solving these two equations for c and A, we have 

A = —3 M 2 (H 2 + H^j , 

c = M 2 a 2 H. (3.5) 

At the first glance, one may worry that temporal diffeomorphism invariance is broken since c and A are 
time dependent. However, the time dependence of the coefficients is not sufficient to break temporal diffs. 
One way to check is to perform a general coordinate transformation x 11 —> x 11 + and then promote all 
four parameters into fields, —> ir^lx). These four 7r M are the would-be Goldstone bosons associated 
with broken diffeomorphisms. It is easy to check that there is no dynamical Goldstone boson associated 
with temporal diffeomorphism breaking, or in other words, temporal diffeomorphisms t —> t + remain 
gauge redundant. 

We turn now to the fluctuation operator. The simplest form of operator that starts linear order in 
fluctuations is 

n ik n kj 

■ (3.6) 

To distinguish it from the metric fluctuation 5g ^ v , we put a bar over 5. The trace vanishes up to linear 
order in fluctuations, b\g u = 0, where the subscript 1 on <5 denotes linear order in the expansion. 

The term quadratic in fluctuations can be obtained equivalently as the second order of the operator 
Ylib 2 g u or the operator bg %3 5g 13 . The two differ only by a factor a which can be absorbed by 
redefinition of the coefficient, and we do not need to write the Yli b 2 d n term. 

The construction of the operator Eq. (3.6) is not unique. For instance, we can equally write 


3^9’v¥+2^9vy- [j>”) ’ - ( 3 - 7 ) 

i,j,k i,j,k \ i J 

All of these operators as well as the products among them, give rise to exactly the same quadratic action 
as dg^bg 13 , up to a prefactor of scale factor a, which can be absorbed into the coefficient. The reason 
is quite simple. Suppose we have a general operator T 13 constructed out of only the spatial metric g 13 . 


E n ik n kl n lm n nr< 

k,l,m9 9 9 9 

n kl n lm n mk 

2-^k,l,m 9 9 9 


g 13 - 3 


Efc 1 g ik g kl g lj 


\ ' n kl n kl 

2^k,l9 9 


g 13 - 3 





Expand the spatial metric around the FRW background and then decompose the linear perturbations 
into trace part and traceless part, 

g ij = a~ 2 5 ij + \8g kk S ij + 5g ij , (3.8) 

o 

where 5g lJ is traceless. The trace part of the perturbation can be absorbed by redefinition of the scale 
factor a. Therefore, if we demand that the background part of operator T 13 vanishes, the trace part of 
the perturbation must also vanish, and we have T 13 = f(a)5g 13 at linear level, where /(a) is a generic 
function of scale factor a. The quadratic order operator can be constructed out of the product of two 
T 13 s, i.e. 6g 13 8g 13 . On the other hand, note that the trace T = T M vanishes at linear perturbation 
level and at non-linear level T ~ 5g l3 5g 13 again. We have thus proven that 5g l3 5g 13 is the only quadratic 
order operator needed. 

Putting these elements together, we write the action Eq.(3.2) as 


5 


d 4 Xy/—g 


-MplZ — 3M 2 


(h 2 + H^j 


+ M 2 a 2 Hg v 


M 2 Mf5g i3 ~5g i3 + ... 


(3.9) 


where the dots stand for the operators starting from cubic order in fluctuations. For the simplicity of 
notation, we drop all summation symbols, all repeated indices should be summed up by default. Noted 
that Mf could be a generic functions of g 11 and thus time dependent. 


3.3 Action for the Goldstone Bosons 


To exhibit the three Goldstone bosons that nonlinearly restore general covariance, we can perform a broken 
spatial diffeomorphism (so called Stuckelberg trick). These three Goldstone bosons are decomposed into 
1 scalar mode and 2 vector modes, in addition to the 2 tensor modes in GR. In the unitary gauge, the 
graviton “eats” the three Goldstone bosons and becomes a massive spin-2 particle with 5 polarizations: 
one helicity-0 mode and two helicity-1 modes in addition to the helicity-2 tensor modes. 

Under the broken spatial diffeomorphism, x l —> x l = x l + £* and t t = t, the metric g lJ transforms 
as 


- i 11 (*M) = ( 3 - 10 > 

and d 4 Xy/—g is invariant under all space time diffeomorphisms. Recalling the construction of the action 
in the unitary gauge, A and c in eq. (3.2) are functions of g u , and their values are fixed by the background 
Einstein equations to cancel out tadpole terms. This requires that under the broken diffeomorphisms, A 
transforms as 


A —)• A + 


dA 

dt 



1 

2 ! 


dg 3 3 \ 1 d 
dt ) dt 


dA 

dt 





(3.11) 


where 5^g is the variation of g u (x) under the broken diffeomorphisms, i.e. the trace of eq. (3.10). Noted 
we have applied chain rule in the Taylor expansion. The coefficient c transforms in the same way as A. 
After the spatial diffeomorphism transformation Eq. (3.10), the action Eq. (3.9) reads 


S = M 2 / d 4 x^g 


(h 2 + H}+ a 2 Hg^ ' x + ^ 


dx^ 


dx 1 ' 


-3 \ 2H+ H 




(3.12) 
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Then we promote the nonlinear parameters £* into scalar fields, £* —> 7r a <5*, and assign to 7r a the transfor¬ 
mation rule 


ir a (x) -> 7r a (x(x)) = v T a {x) - 8f?(x). 


(3.13) 


With this definition, the Goldstone scalars non-linearly recover general covariance and describe the fluc¬ 
tuation around the FRW background. 

The resulting action for the Goldstones is 


5 = M 2 / d 4 x^g 


H 


a \2 


-n - 3 (H 2 + H) + a 2 H (/ + 2 5g a ^ a + g^d^d^) - - ( 2 H + -) ( d a ir a ) 


-2M%d a ir b 5g ab - 


2Mn 


diTr a diir a + ( d a n a ) 2 ) + ... 


(3.14) 


One should keep in mind that the 7rs are the physical excitations over the 3 scalar fields’ vevs in Eq. (2.1). 
Note that in the increased symmetry of the de Sitter limit H —> 0, the kinetic term (7f a ) 2 vanishes and 
the expansion breaks down. 

When the energy scale is much greater than the mass of gauge boson, the mixing between longitudinal 
and transverse components of the gauge field becomes irrelevant (helicity is approximately conserved). 
The two sectors decouple and analysis is greatly simplified in the Goldstone language, where there are 
only interacting scalars. As seen in Eq. (3.14), the leading order mixing is determined by H or and 
this simplification is also achieved when the wavelength is small enough. The leading order mixing terms 
between Goldstone bosons and metric perturbations are 


MlH6g aj djTT a , and M 2 p M%d a K b 5g ab . (3.15) 

We canonically normalize the action for Goldstone pions and metric perturbations as 7r“ ~ M p aH 1 / 2 TT a 
and Sg'c ~ M p 5g i:) , after which the mixing terms read 

H^Sg^djTT* , and . (3.16) 


Note that it is appropriate to compare 3-momentum scales here due to violation of Lorentz invariance. 
We can see that mixing can be neglected for energies above both of H 1 ^ 2 and M|ET _1//2 . As expected, 
the mixing scale is essentially the mass of dynamical modes on the metric perturbations derived below, 
as the co-factor e 1 ' 2 is absorbed in canonical normalization. After neglecting the mixing, the action of 
Goldstone Bosons dramatically simplifies to 


S = M; / d A Xy/^j 


a 2 Hg^d^ a d„K a - M 2^ + f ) (4vr“) 2 - ^ 


diTT a diTT a + ^ {d a TT a ) 2 ) + ••• 


(3.17) 


Away from the short wavelength approximation, for H/H 2 ~ 0(1) (also far away from de Sitter), the 
first coupling term in Eq. (3.16) is important, and to understand dynamics of modes at the momentum 
scale of interest k 2 ~ a 2 H, a full perturbations analysis is necessary, which will be presented in the next 
subsection. In the de Sitter limit H —> 0, the second mixing term diverges, another manifestation of the 
strong-coupling problem. We shall see below in Sec. 3.6 that strong-coupling is avoided with the inclusion 
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of higher derivative terms. If the coefficient M 2 is independent of or only weakly dependent on the scale 
factor, the mixing can become negligible at late time. 

Before moving to the perturbations analysis, we take advantage of the ease with which the breakdown 
scale can be estimated in the Goldstone language. Since Lorentz invariance is broken, A can be written as 
a cutoff in 3-momentum or in energy, which are related approximately by a factor of the sound speed at 
large k. Reading from Eq. (3.14), the interaction terms should be subdominant compared to the kinetic 
energy for the theory to be perturbative. The most stringent constraint on the breakdown scale comes 
from looking at the 3-7T coupling inside the 5g lJ 8g^ term. The 1 —> 2 amplitude becomes of order 1 
at the 3-momentum scale 


A ~a 7 / 2 ^ 3 / 4 


M 


1/2 


M 2 


a 7 / 2 e^ A H 


(HM p ) 1 / 2 

M 2 


(3.18) 


where e = —H/H 2 is the usual slow roll parameter. Provided M 2 is not too close to M p \, this is 
parametrically higher than the mixing scales given in Eq. (3.16). On the other hand, for the effective 
theory to be useful in the cosmological context, the breakdown scale should be (much) larger than H so 
that it effectively describes long wavelength near-horizon dynamics. The requirement A 3 > H is equivalent 
to 


M 2 < a 7 / 2 e 3 / 4 (HM p ) 1 / 2 . 


(3.19) 


Note that other operators, for example Sg^Sg jk 5g kl at third order, also contribute to the 3 - 7 T coupling. 
A priori, the corresponding coefficient M 2 (in obvious notation) is the same order of magnitude as M|. 
However, when expanding the operator to obtain the ti-tt interaction terms, each power of 5g brings with 
it a factor a~ 2 . Consequently, the effect of operators that appear at higher order in unitary gauge is 
typically diluted faster in the expansion of the universe. Only if the associated coefficients scale with 
a compensating power of a can these higher orders be relevant or worse lead to strong coupling at late 
time. As we do not observe a phase transition in the gravitational dynamics (except possibly the end 
of inflation), we may exclude this possibility and consider that terms higher than the M | term are 
suppressed. 


3.4 Full Perturbations Analysis 

We have learned that away from short wave length limit, generally the coupling between Goldstone bosons 
and gravity cannot be omitted. To perform the full perturbations analysis, we first decompose the modes 
according to helicity and then identify and integrate out non-dynamical metric degrees of freedom. This 
procedure results in the effective actions for the 5 dynamical modes, though the calculations are somewhat 
tedious. 

Due to the SO(3) rotational symmetry of our background space time, we can decompose the metric 
perturbations into scalar modes, vector modes, and tensor modes. The helicities completely decouple at 
linear perturbation level. We define the metric perturbations as follows, 


9oo — — (1 + 2(/>) , 


90i 

9ij 


a(t) (St + di/3) , 
a 2 (t) 


5ij + 2 i/j5ij + (didj — -d 2 )E + —(diFj + djFi) + 7 y 


(3.20) 
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where 4>, (3, i/j, E are scalar perturbations, S l7 Fi are vector perturbations and 7 ij is the tensor pertur¬ 
bations. Vector modes satisfy the transverse condition, 

d,S l = d,F l = 0 . (3.21) 

Tensor modes satisfy transverse and traceless condition, 

7 i = da ij = 0 . (3.22) 

Under spatial diffeomorphisms, the vector held defined by [28] 

Z i = ^8 ij {d j E + F j ) (3.23) 

transforms as 

(3.24) 

Comparing to Eq. (3.13), we see that the combination (Z l + it 1 ) is a gauge invariant quantity. In the 
unitary gauge, Z l eats ir l , and survives in the linear perturbation theory. This is in contrast to general 
relativity, where both of E and T) are non-dynamical and can be integrated out. 

In this section, we analyse the metric perturbations in unitary gauge, in which we fix vr* = 0. 


3.4.1 Scalar Modes 


In the scalar sector, <j>, fd and ^ are non-dynamical. After integrating them out, we obtain the quadratic 
action for scalar modes, 



dtd 3 k (lC s E 2 - n s E 2 ^j , 


(3.25) 


where 


/C, = 


-k 4 a 5 H 


4 (k 2 - 3 a 2 H^j ' 


a 3 k A 


U, = 


k A H - HH 36 a A H 2 - 21 a 2 k 2 H + k 


12 H k 2 - 3 a 2 H 


+ 


2 k A M 2 
3 a 


(3.26) 


The scalar mode is ghost free, as long as H is negative. On the other hand, the kinetic term vanishes 
in the limit H —> 0, which implies the strong coupling in this background, as found in the previous 
subsection. The scalar action is canonically normalized by defining the field as 


with the result that 


£ = 


—M 2 k A a 2 H\ ''' 

-?-^ E 

2 k 2 - 6 a 2 H ) 


Sf ] = ^J dtd 3 ka 3 (f 2 - w 2 £ 2 ) , 


(3.27) 


(3.28) 
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where 


*2 = 


8Mn 8k 2 M% 

. t~ t 7 + 


36 a A H b e 


6 ,3 


+ 


3a 2 H A k 2 e ( 77 s + 77 — 77 s + 2e 2 — (77 — 22)e — 2) 


a’ 3a 6 H 2 e ' (A ; 2 + 3 a 2 H 2 e) 2 
H 2 k 4 (77 2 + 277(5 + 5) - (677 + 56)e + 16) 
4 (A ; 2 + 3a 2 H 2 e) 2 

The “slow roll” parameters used here are defined by 

H e 


+ 


2 (, k 2 + 3a 2 H 2 e) 2 
k 6 ( 1 + 2 e — 77 ) 


3a 2 (A : 2 + 3 a 2 H 2 e 


,2 > 




F 2t 


77 = 


tfe’ 


s = 


Hr ]' 


(3.29) 


(3.30) 


In the IR regime A ; 2 <C a 2 H 2 e, the dispersion relation of scalar modes can be perturbatively expanded 
with respect to k, revealing a relativistic dispersion relation 

c 2 k 2 


uj 2 ~ 


+ m s , 


(3.31) 


where 


2 _ e V 2 - ? 7 2 - r] + 77 s 

C.q - ± “T o o “t* 


777 2 = AH z e + 


6 
8 


6 e 


8 M| 
3a A H 2 e' 


(3.32) 


though the speed of sound at momenta k > m s differs from the speed of light. The mass and sound 
speed receive both model independent and model dependent contributions. In general M 2 is a function 
of g u and thus scale factor and time dependent. If Mf oc o 4 , the mass of scalar mode will approach to a 
constant value at late times. 

At short distance, k 2 ^$> a 2 H 2 , the dispersion relation simplifies to 
r, k 2 ( 8 Mo 1 + 2e — 77 \ 8M? 


f 8 Ml 

(jJ~ ~ - • [ --— 

“ - a 2 V 3a A H 2 ( 


+ 


) 8 A/f 2 1 

+ - -H 2 [r] 2 + 2 t /(5 + 5) + 16e 2 - 2{7r, + 24)e + 16] . (3.33) 


At short wavelength, the leading order of mass term in the dispersion relation can also be neglected. 
After doing so, this result agrees with the one calculated in the Goldstone gauge, i.e. eq. (3.17) under 
the helicity decomposition 7 r* = + A 1 where <p is the scalar mode considered just now in the unitary 

gauge. 

3.4.2 Vector Modes 

In the vector sector, Si is non-dynamical. After integrating it out, the effective action for the vector 
degrees of freedom reads 


Si 2) = Ml / fC v F,F l - Q v FiF i , 


(3.34) 


where 


1C V = 


-k 2 a 5 H 


4 (A : 2 - 4 a 2 H^j ’ 


1 2 3 tt M%k 2 

£l v = —~k 2 a 3 H H-1—. 

4 2 a 


(3.35) 
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Similar to the scalar modes, the vector modes are free from ghost instability when H < 0. We define the 
canonical held variable, 


(~M*k*a*H\ 1/2 
\ 2 k 2 - 8a 2 H ) 


and the canonically normalized action is 

= X - J dtd 3 ka 3 (PiF - uj 2 v FF) 


(3.36) 


(3.37) 


where 


wl = 


8M| | ‘2k 2 M% | 64a 4 i7 6 e 3 | 2a 2 H' i k 2 e (r/ 2 + g - rjs + 2e 2 - (77 - 24)e - 2) 

a I a rro I 7“ (7TT7 7"7 — o 


a 6 H‘ 2 e (Aa 2 H 2 e + k 2 f 
H 2 k 4 (rj 2 + 2 r){s + 5) + 8e 2 - 2(5 rj + 36)e + 16) 


+ 


(4 a 2 H 2 e + k 2 Y 
k e 


(3.38) 


4 (4a 2 H 2 e + k 2 ) 2 a 2 {Aa 2 H 2 e + k 2 ) 2 ' 

In the IR regime k 2 <S a 2 H 2 e, the dispersion relation of vector modes is also approximated by the 


relativistic form 


wl - 


elk 2 


+ m 2 v , 


(3.39) 


where 


9 e 77 2 — 71 — 77 2 + 77S 2 Ml 

cl = 1 + - - 7 - '—P -- + 2 


ml = 


8 

8M 9 2 , u2 

—f^+AH 2 e. 
a 4 


8e 


a 4 77 2 e’ 


(3.40) 


Comparing to Eq. (3.32), we see that in the IR limit, the mass of vector modes is the same as that of 
scalar modes. 

At short distance k 2 3> a 2 H 2 , the dispersion relation of vector modes is simplified as 
k 2 


Uv - — 
a A 


2M? \ 8M? 1 9 r 9 / 9 . . n . . 

1 + + [ f/ + 2r?(s + 5) + 8e " 10(??+4)e+ 16 ] • (3 - 41) 


Note that sound speeds of scalar modes and vector modes are not independent. In the short distance 
k 2 3> a 2 H 2 , they are related by c 2 ~ | (1 + c 2 ), under the approximation that all slow roll parameters 
are much smaller than unity, and in this limit agrees with the result of [50] . This is due to the uniqueness 
of the quadratic operator 5g l2 5g 13 at linear perturbation level. 


3.4.3 Tensor Modes 


The quadratic action for tensor modes reads 


o(2) _ M P 


7 ijf J ~ 


+ 


8 M 9 , 


+ 4 He) -fiji 




(3.42) 


The tensor modes also become massive, with the same mass as the scalar and vector modes in the IR 
regime. In the observational aspect, the non-vanishing mass gap leads to a sharp peak on the stochastic 
gravitational waves spectrum. The position and height of the peak carry information on the present value 
of the mass term, as well as the duration of the inflationary stage [57]. 
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3.5 Higher Order Derivatives 

In this subsection, we expand up to second order derivatives. In general, we may expect higher- 
derivative operators to arise at least after calculating loop corrections to the action [65]. To this order, a 
(over) complete set of operators is 

Ai-n, A 2 -K l \ A 3 -U ijij , A 4 • V p V u x a V p V v x b 6 ab , As • Dx a Dx b 6 ab , a 6 -v p a 7 v p a 8 , 
A^V'MioVMn, Aio • V‘ t V'x a V p V ff x b 8 ab ■ f p „f ptT , A n • x a V^V ff x b 5 ab • f pp f„ a , 
Ai 2 f p , u ('V p f pa )('V u x a )'V p 'V a x b 5 ab , A 13 (V p f pa )(V p f p „)(V'x a )(V T x b )6 ab , d^d^, (3.43) 

where A n s are generic functions of g pu and f pv = d p x a d u x b 5 ab . Compared to first order derivative terms, 
these terms are suppressed by the UV scale A~ 2 and thus less relevant at low energy. 

By means of the following metric field redefinition, a generalised conformal transformation, 


9[iv t (1 + B ) g pu — g P vi 


(3.44) 


a theory with second order derivatives is mapped to a theory with first order derivative terms plus third 
and higher derivatives 


d A Xy/^g 


-MplZ + F (< f\ V) + G (V p V"x a ,...) 


d^xyf^g 


l -M 2 p n + F{r v ,f l 


[IV) + 0(V 3 x) 

(3.45) 


where B and G are the functions of higher order derivative terms and 


B = A(g p " 1 f pu )-G, 


(3.46) 


where A (g pv , f pu ) is a function of g pu and , and its form is decided by the parameters in F (g pv , f pv ). 
This field redefinition works equally well if the next leading derivative terms are third-order, with the 
resulting theory containing only first-order derivatives and fourth- and higher-order derivatives. The 
procedure could be repeated to remove derivatives up to a desired finite order: Starting at n > 2, n-order 
derivative operators can be removed in favor of n + 1-order derivative and higher terms. 

As an example, consider an action with two second-derivative terms, 


J V=9 + M 2 m 2 (c 0 + af + c 2 f 2 + c 3 / 3 ) + M. p m (Vv^V [NvX b 5 ab - 


(3.47) 


where / = g pu fn V . It is equivalent to 




^M 2 77 + M 2 m 2 ( c 0 + c 4 f + c 2 f z + c 3 /‘ 


with 


9 [iv 


, (V M V u x b 5 ab - nx a Ox b 5 ab ) 

M p m (2cq + ci/ — c 3 / 3 ) 


9[iv 


(3.48) 


(3.49) 
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We have used the approximation 


V=~9 


{y^x a W^ v x b 5 ab - Ux a Ux b 8 ab ) 
M p m (2c 0 + ci/ - c 3 / 3 ) 


71 - V^gTl, 


(3.50) 


because terms like E\x a E\x b 8 ab -7l are third order derivative terms, so that they are additionally suppressed 
and we can neglect these terms when truncating at the second order derivatives. The effective action for 
Goldstone bosons is derived from eq. (3.9) with H 2 , H and g 13 replaced by the ones induced by Eq. 
(3.49), and then performing the spatial diffeomorphism transformation shown in Sec. 3.3. 

Notably this field redefinition implies that the sound speeds of the scalar and vector modes (at k 2 ^$> 
a 2 H 2 , Eqs. (3.33) and (3.41) respectively) are modified only by small corrections to the cosmological 
parameters H,e,rj,s due to the change in metric. Since the expressions Eqs. (3.33) and (3.41) remain 
valid, the high energy relation 

c ^-^( 1 + c s ) ( 3 - 51 ) 

is preserved even in the presence of second-order derivative terms. Repeating the metric redefinition 
procedure to remove derivatives terms of any finite order, we find that the relation is a robust prediction 
of the effective theory, valid as long as the underlying derivative expansion is valid. 


3.6 de Sitter and Minkowski limits 

The higher derivative terms become important in the de Sitter H —?• 0 and Minkowskian H = H 2 — > 0 
limit. In the limit, the kinetic term from the lowest dimensional operators vanishes, and kinetic terms aris¬ 
ing from higher derivative operators become leading order. This eliminates the strong coupling problem 
in the de Sitter and Minkowski limits. 

For instance, with higher order derivative term d^8g 13 d^Eg 13 taken into account, in Minkowskian limit 
H = H 2 -> 0, we have 


S = J d 4 x^Tg [ — M 2 M 2 8g ij 8g 1J - d 2 M p M 2 d p 8g ij d»8g ij + ... 
= [ d A Xyf^g 


1 / ^ . .' \ 2 r\ 2 i r\ 2 7 . 1^0 J O O k 


2 d 2 M p M 2 ( di^diki + - (<9*7r*) - <9 2 7 t*<9 2 7t* - -didjir j did k n 


—2 MpM '2 • ( diir a diir a + - (S a 7r a ) 2 ) + ... 


(3.52) 


where d 2 = didjdij. d 2 is an 0(1) positive constant since a priori higher derivative terms may be similar 
in size to the leading term and suppressed primarily by the additional powers of k 2 / A 2 . The Goldstone 
action shows clearly how the goldstones obtain non-vanishing kinetic terms directly related to the higher 
derivative terms. After canonical normalization, we see that instead of strong coupling, the corresponding 
scalar and vector modes become massive in the Minkowskian space-time, with masses m 2 ~ d 2 M p M 2 . 

For the tensor modes, the situation is different: the quadratic term M 2 6g^8g lJ still provides a nonvan¬ 
ishing mass, as seen in the H —> 0 limit of Eq. (3.42). The higher derivative terms only make a correction 
to the mass. For instance, if we include the higher order term d^Eg 13 d^dg 13 in the Minkowskian space-time 
limit, the tensor action is 

S D if" J ( X + TO foi' 3 ~ ~ 8M 2 7*j7* J , (3.53) 
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Canonically normalizing 7 ^, the graviton mass is 


~ 8M|(1 


8d 2 M 2 

M p 


(3.54) 


Consequently, the vector and scalar modes have masses m v ,m s ~ V m p/m 2 ■ rriT ■ Considering horizon- 
scale perturbations k ~ H, these modes are relatively heavy and could be integrated out, having 1 /M 2 p 
impact on tensor-mode observables. On the other hand, we would like to informatively mention that the 
sound speed of tensor mode will also be modified, with more of higher derivative terms in eq. (3.43) 
included. In this case, the sound speed of tensor mode deviates from unity by a factor of M 2 /M p . 


To conclude, in this section we have investigated perturbations on the FRW background. We first derived 
the Goldstone action up to quadratic order, which clearly isolated the strong coupling problem as well 
as a possible resolution by the inclusion of higher derivative terms. Seeing that mixing with the metric 
can not in general be neglected, we then performed a full perturbation analysis in the unitary gauge. 
This analysis exhibited a well-behaved massive spin two particle, with 5 polarizations: one scalar mode, 
two vector modes, and two tensor modes. All helicity modes are massive, and the masses shown to be 
identical in the low momentum regime. The dispersion relations of these 5 modes are fully characterized 
by the parameter set {H, e, 7 , s, M|}. 


4 Several Examples 


4.1 The Minimal Model and Next-to-minimal Model 

The simplest example is obtained by setting M 2 = 0 in the general action Eq. (3.9), and only keeping 
the first three terms. In the <\> a language, this theory corresponds to Einstein gravity and three canonical 
massless scalar fields with space like VEV Eq. (2.1) [43], 

S = J d 4 xy/^g - ^rr^g^d^d^dab - , (4.1) 

where A is the bare cosmological constant. The energy density of the spatial condensate scales as p oc a~ 2 
and its equation of state equals to —1/3. In the linear perturbation theory, after canonical normalization, 
scalar, vector, and tensor polarizations of graviton have the same dispersion relations 


2 _ 2 _ 2 _ k2 , 2 ™ 2 
u s u v u t ^ ^ 


(4.2) 


with the same non-vanishing mass. These dispersion relations are identical due to the 50(3) internal 
symmetry of the scalar fields, which has been imposed to ensure the rotational symmetry of the vev 
configuration. For the same reason, the scalars can be re-decomposed into 3 polarizations: two transverse 
modes and one longitudinal mode, 


= S ai (di<p + At ), (4.3) 

where diA 1 = 0. Due to the 50(3) symmetry, we could rotate longitudinal mode a bit “into” transverse 
modes, and on the other hand transverse modes are rotated a bit “into” longitudinal mode, and leave the 
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action invariant. In the unitary gauge, these transverse and longitudinal modes are eaten by graviton, It 
implies the masses of scalar modes and vector modes of graviton should be the same. 

To see how the effective theory operator Zg^Zg 10 in Eq. (3.14) is related to a specific model, we consider 
as an example a general polynomial of the tensor defined in Eq. (3.1), i.e. = d fl x a d u x b Z a b, f = g pv 
For instance, starting from the theory truncated at cubic order, 


5 = M 2 V 


d A X\J—g 


n 


- m 2 (c 0 + Cl f + c 2 f 2 + + c 3 f + d 3 f p f p J" + g 3 f ■ ) 


, (4-4) 


where cq is the bare cosmological constant. In this case, the coefficient M 2 equals to 


M| = m 2 



3d 3 3g 3 \ 

+ ^T + -T ' 
a z a- ) 


(4.5) 


The dispersion relations of the 5 polarizations of graviton can then be read from Eq. (3.26)(3.32), Eq. 
(3.35)(3.40) and Eq. (3.42) directly. 

Generally, given a Lagrangian with the function of — f \J—gF (g pu , f^ u ), we can calculate M| i n this 
way: we first Taylor expand the Lagrangian around background, 


F{g^J„u) 


Fq + 


6F 

JgiF 


Sg^ + 


1 5 2 F 

2 ! 5g^5gP a 


5g^bg pa + .... 


(4.6) 


Then note that at linear perturbation level, Zg 12 equals the trace-less part of metric fluctuation Zg lJ . 
Finally, we decompose the metric fluctuation into trace part and trace-less part, Zg 12 = Zg 12 + ^Zg kk Z l F 
is identified as the coefficient in front of the sum of the trace-less terms (see the appendix A for more 
details). 


4.2 Generalization to Spatially Non-flat Universe 

Up to now, we have analysed the gravity theory with broken spatial diffeomorphisms in a flat FRW 
universe. It is straightforward to generalize it to a non-flat FRW universe. In this case, the internal 
metric of scalar fields configuration is replaced by one which is compatible with the metric on the non-flat 
spatial slice. 

For a spatially non-flat FRW universe, the space time metric can be written as 


ds 2 = — dt 2 + a(t) 2 £lijdx l dx\ 

where Qijdx l dx :] is the metric on a 3-sphere 

0 _ x , KZuZj m x l x m 
ij ~ + 1 - KZ lmX i X ™' 


(4.7) 


(4.8) 


where K = 1 for a closed universe and K = —1 for an open universe. In the unitary gauge, the tensor 
fn U takes the form which compatible with 3-sphere metric, 


Uu = d^ a dA b G ab m = ( 0 , flij). 


(4.9) 


A possible vacuum configuration for scalar fields is 


r 


= X 


GabiZ^') Z a b T 


KZ ac Z hd Z> c Z> d 
1 - KZ cd (j) c (t> d ' 


(4.10) 
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It is easy to check that the above vacuum configuration are indeed on shell and satisfy the equations of 
motion, 


sTy^ v <t> a + g^d^ b d^ c Tl = 0 , 


(4.11) 


where T^. is the affine connection which derived from the inner metric G a b(4>). The generalization of our 
effective field theory approach to a spatially non-flat universe is quite straight forward. Including spatial 
curvature, the effective action can be written as 


S = Alp 


d A Xy/^g 


-n-3 


H 2 + H + a 2 H - K) n ijg lJ - Mi5g lJ 5g lJ + ... 


(4.12) 


where the quadratic order operator is defined by 

SgV = g ik a kj - 3 9 ''‘ 9 ‘^ Qmi ■ (4.13) 

Non-zero spatial curvature is sufficient to ensure the kinetic term is non-degenerate. This suggests that 
another way to cure the strong coupling problem of massive gravity in Minkowski space is to perturb in 
the direction of non-vanishing spatial curvature. 


4.3 A Self-accelerating Universe 

When we apply our massive gravity theory to cosmology, one of most interesting questions is whether or 
not a graviton mass term can accelerate the cosmic expansion. A similar question was studied in Ref. 
[50] , in which they proposed an inflationary model “solid inflation”, with de Sitter-like expansion driven 
by the vacuum energy of the “solid”, that is the spatial condensate vacuum configuration Eq. (2.1). In 
this section, we provide another way to realize a de Sitter phase. 

We work on the static chart of the de Sitter phase, where the metric takes the form 

ds 2 = —(1 — H 2 r 2 )dt 2 -|- \ dr 2 + r 2 d0 2 +r 2 sm8 2 d(j) 2 , (4-14) 

1 — H z r z 

and H is the Hubble constant of de-sitter space-time. The Einstein tensor reads, 

GO = -3H 2 5% . (4.15) 

In terms of spherical coordinate, the 3 scalars can be written as 

4> a = x a , x a = r(smdcos(f>,sm6sm4>,cos9) , (4-16) 

and the tensor 

fiw = d^d^dab = (0,1, r 2 , r 2 sin(9 2 ) . (4.17) 

Normally, such a field configuration is not compatible with the space time metric Eq. (4.14), since the 
scalars do not satisfy their equation of motion. For instance, if we consider a canonical Lagrangian, 

M 2 J d^xyf^g - ^mV^0“3 1/ /<5 a6 j , (4.18) 
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the energy momentum tensor reads, 



= -H 2 r 2 - 
2 

= -~H 2 r 2 
2 

= - H 2 r 2 - 
2 

_ rri2 

~ 1 2* 


3 

2’ 



1 

2 ’ 


(4.19) 


Comparing to Eq.(4.15), we can see Einstein equations are not satisfied. On the other hand, the vacuum 
configuration does not satisfy the equations of motion for scalar fields 4> a either, 

gT'V lx V v 4> a oc H 2 r ± 0. (4.20) 


The SO (3) scalar fields’ configuration in the static chart of de Sitter phase implies that there are large 
shears and energy momentum flows in the non-static coordinates. Nevertheless, if we include the higher 
order kinetic interaction terms, by tuning the parameters we may be able to cancel out the shears and 
flows and realize a self-consistent de Sitter solution. We found this solution at least contains 4th order 
kinetic interactions, 


M 2 m 2 J d 4 xy/=g [a f + c 2 if 2 + c 22 fUf" + c 3 i/ 3 + c 32 fpf”fP + c 33 f ■ fPf” 

C4l/ 4 + C42 + c A3 f • + cuitffp 2 + C 45 f 2 • /£/; ] , (4.21) 


where all coefficients ci, C 21 , C 22 --- are coordinate independent constants. As a self-consistent solution, 
it must satisfy Einstein equations, T„ = G„ = —3 H 2 5„, and the equation of motion for scalar fields as 
well. These conditions lead to a set of nontrivial constraints on the coefficients, 


C21 

C22 

C41 

C42 

C43 

C44 

C32 

C31 


- (—3ci - 2 c 33 ) , 

3 

- (ci - 2 c 33 ) , 

ig c i - h ^+ 4c45 ) ’ 

3ci c 33 

16 +C45 8 ’ 

1 , , 1 

-(C33-8C45)-—ci, 

- (-Cl - 4 c 4 5 ) , 

-c 33 , 

0, 


(4.22) 


and the energy momentum tensor can be calculated as 

Ty = ^ (3ci - 2 c 33 ) M 2 m 2 5t . (4.23) 

Indeed with this choice of parameters, the vacuum energy behaves as an effective cosmological constant. 
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However, as we mentioned in the section III. B, in the de Sitter limit, Goldstone bosons become 
strongly coupled, and the theory breaks down. To show how it occurs in this model, we perturb the 
scalar fields’ configuration by introducing the (Goldstone) excitations 

<j> a = x a + ir a . (4.24) 

Then we take the decoupling limit to decouple Goldstone bosons and gravity, 

m —¥ 0, M p —» oo, A 2 = M p m = const. (4.25) 

In momentum space, the quadratic action for Goldstone bosons reads 

4- 2) = J ^ ab n a n b - c 2 diag k 2 TT a ir b 5 ab - c 2 mix (k a ir a ) 2 , ( 4 . 26 ) 

where 

/C = — (ci + 6 c 2 i + 2c 22 T 27C31 + 3C3 2 + 9C33 + IO8C41 + 4C4 2 + I2C43 T I2C44 + 36C45) , 

Cdiag = — (c 1 + 6c 2 i + 4c 2 2 + 27 c 3 i + 9C3 2 + 15 C 33 + IO8C41 + 16 C 4 2 + 3OC43 + 24 C 44 + 54c4s) , 

C mix = ^ (2 c 2 i + C 22 + I8C31 + 3C3 2 + 7C33 + IO8C41 + 6042 + 15 C 43 T 14 C 44 + 39C4s) . ( 4 . 27 ) 

After inputting the constraints on these coefficients, i.e. Eq. (4.22), we find that JC = c diag = = 0 

and the quadratic action vanishes identically. 

As we saw above in the effective theory analysis, this degeneracy arises in the exact H = 0 limit, and 
one of possible solutions to the strong coupling problem is to introduce a small deviation from de Sitter 
space time. Another possible solution is to include higher order derivative terms. As discussed in section 
3.6, the higher order derivative terms give rise to a non-vanishing kinetic term for Goldstone bosons even 
in de Sitter spacetime and thus heals the strong coupling problem. 


5 Conclusion and Discussion 


In this paper, we characterize the most general theory of spatial diffeomorphisms breaking. By means 
of effective field theory approach, after writing down all possible operators compatible with underlying 
symmetries, we are able to describe the effective theory of fluctuations around the FRW background with 
spatial diffeomorphisms breaking. We showed that the most generic action on a FRW background can be 
written in the form 


S 




-M 2 n - 3M 2 (H 2 + H 


+ M 2 a 2 Hg 11 


M 2 Mf5g^S f f + ... 


(5.1) 


where 5g lJ is the covariant operator constructed out of g l] and starts from linear order in fluctuations. 
Differences among models (UV completions of the theory) are encoded in this operator as well as higher 
order operators. The three broken spatial diffeomorphisms acquire three Goldstone bosons. The cou¬ 
plings between different helicity modes are characterised by the scales H and Ad^/H. Above the mixing 
scale, we can neglect the couplings between different helicity modes, and the action of Goldstone bosons 
dramatically simplifies to 


S = M 2 / cfxyf^g 


1 


H 


a 2 Hg^d^ a d v v a - - \2H + - (d a ir a Y - 


. 2 2M? 


diir a diTr a + ^ {d a n a ) 2 ) + ... 


(5.2) 
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In the unitary gauge, the 3 Goldstone bosons are eaten by graviton, and the graviton becomes a massive 
spin-2 particle with 5 well-behaved degrees of freedom that are decomposed into 1 scalar mode, 2 vector 
modes, and 2 tensor modes. We performed the cosmological perturbation calculation, derived the effective 
quadratic action for each of the 5 polarizations, and found that at linear perturbation level, all 5 polar¬ 
isations have non-vanishing masses. The dynamical properties of these 5 polarisations are characterised 
by the parameter set {H, e, 77 , s, M 2 }, where e, 77 , s are “slow roll” parameters. 

With only first-order derivative operators, the kinetic terms of these three Goldstone bosons vanish 
in the de Sitter, as well as Minkowski limit, where H —> 0. In this limit, the kinetic terms arising from 
higher derivative operators become leading order. In this case, the three Goldstone bosons appear to be 
supermassive at the low energy scale, and thus can be integrated out, leaving two massive tensor modes 
as the only propagating degrees of freedom. 

Away from de Sitter and Minkowski limit, in the IR regime where k 2 <C a 2 H 2 e, the quadratic action 
of the 5 polarisations in the unitary gauge reads 

<S (2) = ^ J dt,d 3 ka 3 t 2 - + m g^ £ 2 + FiF ~ + 7u'7 U “ {^2 + m l) 7ij7 u , 

(5.3) 

where £ is the scalar mode, J- t is the vector mode, is the tensor mode. The graviton mass m 2 = 

4 H 2 e + takes the same value for all 5 polarisations in the long wavelength limit. The sound speeds 

of scalar and vector modes are functions of “slow roll” parameters. 

Several examples on the applications of our formalism to cosmology are also presented. We started 
from the example of a flat FRW universe and a generic polynomial of derivative couplings and showed how 
to determine the dispersion relations of the 5 graviton polarizations. We also generalized to a spatially 
non-flat FRW universe and a static chart of de Sitter space time. 

Along this line, there are several possible extensions of our formalism. For instance, it would be nec¬ 
essary and more realistic to consider the matter distribution in the universe and develop a new effective 
theory with the coupling to mattertaken into account. On the other hand, if we go beyond linear pertur¬ 
bation theory, the interactions among three pions can be characterized by introducing the higher order 
operators like 5g i] 5g 3k 8g kl . This type of non-trivial interaction leads to non-Gaussianity in the metric 
fluctuations. It would be very interesting to study its possible imprint on CMB and large scale structure. 
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A Calculating Mf for a given Lagrangian 


In the section 4.1, we pointed out that for a generic Lagrangian with spatial diffeomorphisms breaking in 
the unitary gauge, 

s D -Ml J y/=ijF 

= -Ml J {f) , (A.l) 

the associated quadratic operator in the effective field theory can be calculated by tracking the traceless 
part of metric perturbation. To do so, we decompose the metric as in Eq. (3.8), 

g ij = a~ 2 6 ij + \dg kk 5 ij + 5g ij , (A.2) 

o 

where 5g 13 is the traceless part of metric perturbations. We expand the action Eq. (A.l) perturbatively 
in terms of the above metric decomposition, 

5 D -«? / ^ (f<«> + i + -) 

d Ml J F a r F m SgW i . + ... 

= Ml jH%‘r Sg ii _ + - (A.3) 

where we have used the 0 — 0 component of Einstein equation 3 H 2 = F in the last line of the above 
formula. 

On the other hand, the effective action for the traceless metric perturbation in terms of our effective 
field approach reads 


S D Ml / -3 [H 2 + H) + a 2 Hg u - M%5g ij 5g ij + ... 


D M 2 / --a'H z 5g l3 5g 13 - a 6 M^bg l °5g lJ + ... 


(A-4) 


By comparing eq. (A.3) (A.4), we find 


M 2 m 9 ij = w* 


2 SgvSg 


(A.5) 
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